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We consider the self-similar solutions associated with the 
critical behavior observed in the gravitational collapse of 
spherically symmetric perfect fluids with equation of state 
p = a/i. We identify for the first time the global nature of 
these solutions and show that it is sensitive to the value of a. 
In particular, for a > 0.28, we show that the critical solution 
is associated with a new class of asymptotically Minkowski 
self-similar spacetimes. We discuss some of the implications 
of this for critical phenomena. 



I. INTRODUCTION 

One of the most exciting developments in general rel- 
ativity in recent years has been the discovery of criti- 
cal phenomena in gravitational collapse. For a variety 
of spherically symmetric imploding matter fields, there 
is a self-similar critical solution containing a naked sin- 
gularity which separates models which collapse to black 
holes from those which disperse Sometimes a dis- 

crete similarity is involved but, in other circumstances, 
the critical solution seems to be represented by a con- 
tinuously self-similar model. This is one which has a 
homothctic Killing vector and contains no dimensional 
constants. Perfect-fluid models of this kind necessarily 
have an equation of state of the form p = afi and so only 
in this case could the critical solution be homothetic. 

Self-similar spherically symmetric perfect-fluid solu- 
tions have been much studied in general relativity (see Q 
and references therein) and the attempt to understand 
critical phenomena has led to a several further studies 
^-^|. However, their precise relationship with the critical 
solution has remained obscure. This is mainly because 
the full family of such solutions had not been identified 
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when critical phenomena were first discovered. However, 
recently Carr & Coley [|o) have presented a complete 
asymptotic classification of such solutions. Furthermore, 
by reformulating the field equations for these models in 
terms of dynamical systems, Goliath et al. |ll| , p^| have 
obtained a compact three-dimensional state space repre- 
sentation of the solutions and this leads to another com- 
plete picture of the solution space. These investigations 
have resulted in the discovery of a new class of 'asymp- 
totically Minkowski' self-similar spacetimes. 

In this paper we shall discuss these new solutions and 
show why they are intimately related to critical phenom- 
ena ||. We thereby demonstrate for the first time the 
global nature of the critical solution. Although the de- 
tailed derivation of these new solutions is given elsewhere, 
this is the first published announcement of their existence 
and the first attempt to link them to critical phenomena. 
The purpose of this paper is therefore to highlight this 
result in advance of the more extensive analyses (since 
these are in a much broader context). In particular, we 
will show how the global features relates to the equation 
of state parameter a and explain why there is only one 
such solution for each a. It should be emphasized that 
numerical studies of the critical solution are always re- 
stricted to some finite range of the self-similar variable z. 
However, as the critical index is approached, the extent 
of the self-similar region grows and one could in principle 
go to arbitrarily large values of z. This paper can there- 
fore be regarded as predicting the characteristics of these 
solutions. 

The discussion will mainly be in terms of the com- 
pact state space |ll],[l2| but, to extract important physi- 
cal features, some of the quantities used by Carr & Coley 
JlOt will be plotted. No equations will be used because 
the discussion is intended to be purely qualitative and 
thereby accessible to the general reader. However, some 
technical terms will be used in the next section, so we 
here give some background references. For an introduc- 
tion to dynamical systems theory in general relativity, see 
|l3| ; for its particular application in the spherically sym- 
metric context, see jn], 12[; for more details of the other 
types of self-similar solutions, see Carr et al. Qj. At the 
end of the paper, we emphasize our key predictions, so 
that "critical" workers can investigate these. 



1 




FIG. 1. The three-dimensional state space of self-similar 
spherically symmetric perfect-fluid models for a > 1/5, ob- 
tained by matching a spatially self-similar region with two 
timelike self-similar regions along the lines H [juj. Labelled 
equilibrium points act as asymptotic states for orbits. The 
triangular surfaces indicated by dashed lines are sonic sur- 
faces. Two of these surfaces contain a sonic line, indicated by 
a line of equilibrium points. 

II. THE SOLUTION SPACE 

We shall focus on spherically symmetric self-similar 
solutions in which the spacetime admits a homothetic 
Killing vector. This means that all dimensionless vari- 
ables depend only on the self-similar variable z = r/t, 
where r is the comoving radial coordinate and t is the 
time coordinate. In a dynamical systems approach these 
solutions correspond to orbits in a three-dimensional 
compact state space. The state space of self-similar 
spherically symmetric perfect- fluid models (for a > 1/5) 
is presented in Fig. [IJ A point in this space corresponds 
to a certain geometry and matter field configuration on 
a homothetic (constant z) slice, while an orbit in the 
state space represents an entire spacetime. All continu- 
ous orbits are future and past asymptotic to one of a few 
solutions with higher symmetry. These appear as equi- 
librium points on the boundary of the state space and 
these points are labelled in Fig. [j]. A physical descrip- 
tion of the solutions asymptotic to them is given in Table 

B 

The state space is divided into two halves, one cor- 
responding to positive z, the other to negative z. This 
means that the solutions in one half are the time-reverse 
of solutions in the other half, so all equilibrium points ap- 
pear twice in Fig. |l|. The sonic surfaces are also depicted 
and solutions generally develop a shock- wave here [ p"5| |. 
However, in two of the sonic surfaces there is a sonic 
line and solutions which pass through this line can be 
extended continuously through the sonic surface. Only 
these solutions will be considered to be physical and the 
number of such solutions is strongly restricted. 

The state space has the advantage that it gives a pic- 
torial representation of the relationship between different 
solutions and the connection between the initial and fi- 
nal states, thereby yielding insights into the global nature 



Label Interpretation 

C Solutions with a regular center or 
infinitely dispersed solutions 

M, M Asymptotically Minkowski solutions 

K Non-isotropic singularity solutions 

F Asymptotically Friedmann solutions 

T Exact static solution 

TABLE I. Interpretation of solutions asymptotic to the 
given equilibrium points. 

of the solutions. However, it has the disadvantage that 
it is rather abstract. To better understand the physical 
aspects of the solutions, it is useful to consider some of 
the physically interesting quantities which arise in the co- 
moving approach. The dependence of these quantities on 
z corresponds to two-dimensional projections of orbits in 
the full state space. Following [jl(|, we use: (1) the scale 
factor S, which fixes the relation between the comoving 
radial coordinate r and the Schwarzschild radial coordi- 
nate R — Sr, this indicating when a solution expands 
infinitely (5 — > oo) or encounters a singularity (S — > 0) 
for finite values of z; (2) the velocity V of the spheres 
of constant z relative to the fluid, which is important for 
the identification of event horizons (\V\ = 1) and naked 
singularities, see [|6|; (3) the density profile fj,t 2 , which 
gives the matter distribution at a given comoving time t\ 
and (4) the mass function 2m/ R, where m(r) is the mass 
within radial coordinate r, this indicating the presence 
of an apparent horizon (2m / R =1), see |p7| . 

We shall first briefly review the previously known fam- 
ilies of solutions. All of these solutions are discussed in 
more detail in where the dependence of the above 
functions on z are shown explicitly. We shall then con- 
sider the new asymptotically Minkowski solutions. This 
is the first discussion of their physical features and the 
first detailed analysis of their relevance to critical phe- 
nomena. 

Asymptotically Friedmann solutions 
There are two one-parameter sets of solutions that are 
asymptotic to the flat Friedmann solution, all of which 
are connected with one of the Friedmann points F. One 
has positive z and the other has negative z. Two 
qualitatively different families can be distinguished: (1) 
expanding-recollapsing solutions (F-K orbits); and (2) 
ever-expanding (or ever-contracting) solutions (F-C or- 
bits), where C is to be interpreted as an infinitely dis- 
persed state. The latter family contains the flat Fried- 
mann solution itself. Thus the flat Friedmann solution 
appears both as an equilibrium point F and as an orbit 
in state space, corresponding to different slicings. 

Asymptotically quasi-static solutions 
For each value of a, there is a unique static solution, 
originally found by Tolman [fL8f . The corresponding (T- 
T) orbit traverses the entire state space and spans both 



2 



positive and negative z. Furthermore, there is a two- 
parameter set of solutions with behavior resembling the 
static solution at large \z\ (i.e. at early times). They 
are all associated with K points, corresponding to non- 
isotropic singularities. As with the asymptotically Fried- 
mann solutions, there are two different families within 
this class: (1) expanding-recollapsing solutions (K-K or- 
bits); and (2) ever-expanding (or ever-contracting) so- 
lutions (K-C orbits). The latter contain the naked- 
singularity solutions discussed by Ori & Piran |l6f| and 
Foglizzo & Henriksen |p|. Unlike the asymptotically 
Fricdmann solutions, the asymptotically quasi-static so- 
lutions necessarily span both positive and negative z. 

Asymptotically Minkowski solutions 
When a > 1/5, solutions exist that are 'asymptoti- 
cally Minkowski', in the sense that the state-space or- 
bits asymptote to equilibrium points that correspond 
to Minkowski space. There are actually two subclasses 
of such solutions, associated with different equilibrium 
points in state space: Class A solutions are connected 
with the M points and are described by two parameters. 
Class B solutions are connected with the M points and 
are described by one parameter. Both of these subclasses 
contain two different families of solutions: (1) singular so- 
lutions (K-M, K-M orbits); and (2) regular solutions (C- 
M, C-M orbits). Only the latter contain a sonic point. 
All these types of solutions are illustrated in Figs. | and 
^, where the arrows indicate whether solutions are future 
asymptotic to M or M. 

Solutions in class A have V — > 1, S — > oo, fit 2 — > 
and 2m /R — > at some finite value z = £*. Although 
this limit is reached at finite z, it should be pointed out 
that most investigations of the critical solution use the 
physical distance R, which is infinite. Examples of such 
solutions are illustrated by the dotted curves in Figs. || 
and |[ Solutions in class B have V — ► V* > 1, S — > oo, 
fit 2 —> and 2m /R — > as z — > oo. Examples of these 
solutions are represented by the dashed curves in Figs. 
U and |^. Both classes are asymptotically dispersive and 
solutions asymptotic to K points also have S — > at a 
finite value of z, indicating the formation of a singularity 
in the past (assuming the time direction indicated in the 
figures). Even though 2m /R — » for both classes, the 
mass m need not vanish. 

It should be emphasized that these solutions are not 
asymptotically flat in the usual global sense, in which 
there is a certain radial decay of the curvature towards 
spatial infinity, see, e.g., pfj| ]. We are not considering 
an isolated system here but rather a fluid spacetime in 
which the Minkowski geometry is obtained asymptoti- 
cally along certain coordinate lines. It can be shown 
that the curvature vanishes asymptotically as the M (M) 
point is approached and r — > oo (and hence t — > oo). Be- 
cause the fluid becomes infinitely diluted, the situation 
is analogous to that of the open Friedmann solution, in 
which the Milne solution is approached asymptotically 




FIG. 2. Asymptotically Minkowski solutions of class A 
(dotted) and class B (dashed) as orbits in state space. Ar- 
rows go from C or K towards the infinitely dispersed state. 
Densely dotted and short-dashed curves correspond to reg- 
ular solutions, while sparsely dotted and long-dashed curves 
correspond to singular solutions. 
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FIG. 3. Physical quantities for the asymptotically Min- 
kowski solutions. The dash-dotted line in the V(z) diagram 
corresponds to the sonic surface. Other designations are given 
in the caption of Fig. B. 



along certain time- lines at late times (see, e.g., Pjl). 



III. THE CRITICAL SOLUTION 

Critical phenomena in gravitational collapse were first 
studied by Choptuik |3] and remain an active field of 
research, see, e.g., J2lj and references therein. The so- 
lution at the threshold of black-hole formation in spher- 
ically symmetric radiation fluid collapse, corresponding 
to a = |, was studied by Evans & Coleman J5J. They 
found it to be a self-similar solution distinguished by the 
following criteria: 

(1) It is everywhere analytic, or at least C°°. In par- 
ticular it has a regular center, and also crosses the sonic 
surface in an analytic way. 

(2) It has a collapsing interior surrounded by an ex- 
panding exterior. This means that the radial fluid three- 



3 



velocity Vr associated with a Schwarzschild foliation 
(which is different from the function V) has exactly one 
zero. 

Subsequently, other authors have used these cri- 
teria to investigate the critical solution for other values 
of a. For a recent review, see [f2~l| . 

The uniqueness of the critical solution can be under- 
stood as follows: For each value of the equation-of-state 
parameter a, there exists a one-parameter set of solu- 
tions with a regular center and a one-parameter set of 
solutions analytic at the sonic line ( Q). Thus it is not 
surprising that the first condition leads to a discrete set 
of solutions. There is only one solution in this set that 
satisfies the second condition and this is the critical so- 
lution. 

We now examine the critical solution in terms of both 
the state space of the self-similar spherically symmetric 
perfect-fluid solutions and the behavior of the various 
physical quantities. The results are summarized in Figs. 
| and | 

Starting from the regular center C, a numerical inves- 
tigation shows that for all equations of state, the orbit 
of the critical solution passes through the sonic line and 
enters the spatially self-similar region (with \V\ > 1). It 
turns out that for a in the range < a < 0.28, it is of the 
asymptotically quasi-static kind: it passes through the 
spatially self-similar region and enters a second timclikc 
self-similar region, finally reaching another sonic point 
(indicated by 'x' in the figures) , which is generally irreg- 
ular. However, this does not invalidate the solution as 
being the critical one, since the solution describing the 
inner collapsing region is usually matched to an asymp- 
totically flat exterior geometry sufficiently far from the 
center. An example of a solution belonging to this class 
is given by the full curves in Figs. ^ and ^. 

We find that for the limiting case a « 0.28, the crit- 
ical solution is an asymptotically Minkowski solution of 
class B, whose orbit ends at an equilibrium point M. In 
Figs. || and |5[ this solution corresponds to the dashed 
lines. For 0.28 < a < 1, we find that the critical solu- 
tion belongs to the asymptotically Minkowski solutions 
of class A, whose orbit ends at an equilibrium point M. A 
solution representing this class is indicated by the dotted 
curves m Figs. |and|. 

For a > a* « 0.89, the investigations in (T^J]]|| indi- 
cate that the critical solution is already irregular at the 
first sonic point. As the matching must be performed 
outside the sonic point [|l6Ul9| , the solution would then 
be unphysical. However, Neilsen & Choptuik [^2| have 
recently demonstrated the existence of a regular critical 
solution for a > a» as well. Our present investigation 
supports their analysis. 

To understand what happens for a = a*, we consider 
equations of state near this value. It turns out that the 
behavior can be understood in terms of the stability near 
the sonic line. In order for a solution to be regular at the 
sonic surface, the corresponding orbit must approach the 
sonic line along one of (at most) two possible directions. 




FIG. 4. Critical orbits in the state space. Arrows go from 
the regular center towards the infinitely dispersed state. Note 
that the orbits correspond to different equations of state, i.e., 
different values of a, and thus belong to different state spaces. 
The three orbits exemplify the cases < a < 0.28 (full curve), 
a x 0.28 (dashed), and 0.28 < a < 1 (dotted). 




FIG. 5. Physical quantities for the critical solution. The 
dash-dotted lines in the V(z) diagram correspond to the sonic 
surfaces. Other designations are given in the caption of Fig. 

Note that for a > 0.28 (dotted), the critical solution tends 
to a finite value of z. 
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Each of these directions is associated with an eigenvalue 
- the direction corresponding to the smaller eigenvalue is 
called dominant and is associated with a one-parameter 
family of solutions (containing just one C°° solution), 
the other is called secondary and is associated with an 
isolated solution. For a < a*, the critical solution cor- 
responds to the secondary direction. However, when a 
gets close to a* , the eigenvalue associated with the crit- 
ical solution approaches that of the other direction. For 
a = a*, the eigenvalues (and directions) are equal, cor- 
responding to a degenerate node, and for a > a*, the 
critical solution is associated with the dominant direc- 
tion. Thus a transition from the secondary direction to 
the dominant direction occurs at a = a*. This transition 
results in severe numerical difficulties, so very high nu- 
merical precision is required to investigate such solutions, 
as pointed out in |22f|. 



IV. CONCLUSIONS 

Our key predictions can be summarized as follows: 

(1) When one gets sufficiently close to the critical so- 
lution that the large- z behavior can be studied, then this 
solution should have the various asymptotic features we 
predict for different values of a. 

(2) There should be a sudden transition in the nature 
of the critical solution as a passes through 0.28, with the 
solution going from the asymptotically quasi-static form 
to asymptotically Minkowski form. However, it should 
be emphasized that the solution is only flat towards null 
infinity for a > 1/3, so one still needs to match to a 
non-self-similar region on a spacelike surface. 

(3) Although the asymptotically flat limit is reached at 
finite z, most critical workers use the physical distance, 
which is infinite. However, it should be pointed out that 
in the stiff case (a — 1), the asymptotically flat state is 
reached at finite physical distance, which should lead to 
some anomalies, see frL0[ ). 

(4) Although we have not explained why the critical 
solution is analytic at the sonic point (this presumably 
relates to the usual stability criterion), we have used the 
asymptotic features to explain why the analytic solution 
is unique for given a. The relationship to solutions which 
are regular but not analytic at the sonic point is discussed 
in more detail by Carr & Henriksen p3[. 
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